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ABSTRACT 

The paper presents a comparison of emission altitudes in very young and very old 
radio pulsars. The author confirmes that altitudes at which radio emission at a given 
frequency is generated depend on pulsar period and on pulsar age, although the latter 
dependence is quite weak. 
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1 INTRODUCTION 

The radial location of the radio emission regions in pulsar 
magnetospheres was discussed by several authors. Using var- 
ious methods reviewed by Cordes (1992, 1993), constraints 
on the emission radii were obtained in a number of papers. 
These estimates show that the pulsar emission region is rel- 
atively compact and lies near the surface of the neutron 
star, apparently closer than about 10 per cent of the light- 
cylinder radius (Cordes 1978; Blaskiewicz, Cordes & Wasser- 
man 1991; Phillips 1992; Gil & Kijak 1993; Hoensbroech & 
Xilouris 1997; Kijak & Gil 1997, 1998; Kramer et al. 1997; 
Gupta, Bhat & Rao 1999). 

Assuming that pulsar is a rotating magnetic dipole we 
can calculate its characteristic age r = P/2P (Manchester 
& Taylor 1977), which is a good estimate of the true age 
providing that P > Po, where P and Po are current and 
initial spin periods, respectively. Most pulsars are 10 6 - 10 7 
years old, however, one separates young and old objects with 
characteristic ages r about 10 4 yr and 10 9 yr, respectively. It 
is important to note that, the group of old pulsars contains 
not only millisecond pulsars, but also pulsars with typical 
period P ~ 0.5 s (e.g. PSR B1848+04, B1952+29) as well 
as recently discovered very old pulsar J2144— 3933 with the 
longest measured preiod 8.5 seconds (Young, Manchester, 
Johnston 1999). It is commonly believed that pulsars form 
two dinstinct populations: high-magnetic field objects that 
are radio pulsars for relatively short times and recycled ob- 
jects which can be very old. 

The millisecond pulsars (MSPs) have different evolu- 
tionary history compared to normal pulsars. They are re- 
cycled pulsars, i.e. they acquired their rapid rotation rates 
(~ a few ms) due to accretion-driven spin-up (Bhattacharya 
1996). Consequently, their magnetic dipole field strength at 
the surface B d = 3.2 x 10 19 (PP) 1/2 G (Manchester & Taylor 
1977) is relatively low (~ 10 8 — 10 9 G), as compared with the 
whole population of radio pulsars (~ 10 12 G). Previously, it 



was thought that MSPs also have different characteristics 
of radio emission from those of normal pulsars. However, it 
has been recently argued (Kramer et al. 1998, 1999; Gil & 
Sendyk 2000) that MSPs emission properties do not differ 
from typical pulsars. 

The semi-empirical relationship proposed by Kijak & 
Gil (1997) describe the dependence of emission altitude on 
pulsar period P (in seconds) , observing frequency v (in GHz) 
and characteristic age r (in 10 6 yr) 



/rr , r iD -0.21±0.07 -0.07±0.03 D 0.33±0.05 

r K G = (55 ± 5)Rv GHz r 6 P 



(1) 



where 7? = 10 cm is the neutron star radius. This was 
obtained from analysis of emission altitudes for a number 
of pulsars at two frequencies. Emission altitudes were esti- 
mated using the pulse width method (see Section 2) based 
on a number of straightforward assumptions: (i) the pulsar 
radiation is narrow-band with a radius-to-frequency map- 
ping operating in the emission region, (ii) pulsar emission is 
relativistically beamed tangentially to dipolar magnetic field 
lines and (iii) the extreme profile wings originate at or near 
the last open field lines. This, of course, involves an implicit 
assumption that the magnetic field in the radio emission re- 
gion is purely dipolar (see Section 3 for discussion). The un- 
certainties given for various quantities in equation (1) follow 
from systematic errors in the pulse width W measurements 
and random errors in estimations of the inclination angle a 
between the rotation and magnetic axes and the impact an- 
gle f5 of the closest approach of the observer to the magnetic 
axis. 

First suggestion that radio emission altitudes at a given 
frequency depend on the pulsar period (r cm oc p~ - 5 ) wa s 
published by Gil & Kijak (1993). This period dependence 
was further carefuly investigated by Kijak & Gil (1997, 1998; 
hereafter Papers I and II) and more precise relationship 
r e m(P) oc p~ - 33 W as obtained. Anyway, the emission re- 
gion is apparently located close to the neutron star in short- 
period pulsars and further away in longer-period pulsars. 
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Figure 1. Emission altitudes at 1.4 GHz versus pulsar period for young (circles) and old (squars) pulsars (data taken from Table 2 
and for old PSRs: B1855+09, J1022+1001 and J1713+0747 from Paper II). Formal weighted fit for pulsars with r 6 < 0.1 gives r cm = 
R(55 ± 8)P°- 25±(U,S and for r 6 > 1000 gives r om = R(35 ± 6)P a3s±a04 . 



Alternatively, the ratio rxG /^lc decreases with increasing 
period (where tlc is the light cylinder radius), what was 
quantitatively demonstrated in Papers I (Fig. 3b) and II 
(Fig. 4a and 4b). The equation (1) describes also a radius- 
to- frequency mapping (RFM) r(v) oc v~ v (where p ~ 0.21), 
which means that different frequencies are generated at dif- 
ferent altitudes above the magnetic polar cap. The RFM 
models of pulsar emission region provide an attractive ex- 
planation of why profile widths (usually) increase with de- 
creasing frequency. Kijak & Gil (1998) studied the concept 
of RFM using the multifrequency data of low-intensity pulse 
widths for normal and millisecond pulsars. They obtained 
RFM with the average exponent p — 0.26 ± 0.09 (for 16 
pulsars), consistent with previous estimates using different 
methods (see Table 4 in Paper II). They ultimately con- 
firmed that the emission altitude depends on frequency v, 
pulsar period P as well as on the pulsar age r, although the 
latter dependence was weak. In this paper we investigate 
a putative weak dependence of emission altitudes on the 
characteristic age indicated in earlier work (equation 1), by 



selecting a sample of very young and very old radio pulsars 
and applying to them the same analysis methods as used in 
Papers I and II. 



2 PULSAR CHARACTERISTIC AGE AND 
EMISSION ALTITUDE 

In this section we calculate emission altitudes for very young 
and very old pulsars using 1.4 GHz data from Effelsberg 
Observatory. In addition, data for two southern pulsars 
B1727-33 (Gould & Lyne 1998) and B0833-45 (Johnston, 
Nicastro & Koribalski 1998) were included because of their 
young age. All these pulsars (Table 1 and Fig. 1) were not 
analysed in previous work (Papers I and II). For estima- 
tions of the emission altitude we use the geometrical method 
(pulse profile width W measurements at low-intensity level 
of the maximum intensity) described in Papers I and II. 

Assuming that the magnetic field in the emission region 
has a dipolar form, the emission altitude r cm (at which radi- 
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W 

Table 1. Emission altitudes for old and young pulsars at 1.4 GHz. We calculated r cm from Eq. 2 and tkg from Eq. 1. Pulse width 
measurements Wq at low-intensity level (i.e. Kb 0.1 per cent) and calculations of opening angle po are presented. References are marked 
for a and (i angles (see Paper I). 



PSR r 6 P r om r KG W p Rcf. 

(10 6 yr) (sec) (r/R) (r/R) (deg) (deg) 



B0833-45 


0.011 


0.089 


30±3 


32 


39.0±1.4 


22.8 ± 1.3 


R93 


B1800-21 


0.016 


0.133 


39±11 


35 


137.±8.0 


21.2 ±4.1 


G94 


B1823-13 


0.022 


0.101 


30±7 


31 


120.±4.0 


21.3 ±3.6 


G94 


B1727-33 


0.025 


0.139 


31±6 


34 


86.0±7.0 


18.5 ±2.7 


G94 


B2334+61 


0.040 


0.495 


59±7 


51 


34.8±2.6 


13.5 ± 1.3 


G94 


B1916+14 


0.089 


1.180 


50±10 


68 


16.2±2.5 


8.0 ±1.2 


G94 


B1848+04 


2818 


0.284 


21±5 


19 


107.±9.0 


10.7 ±2.1 


R93 


B1952+29 


4168 


0.426 


26±6 


17 


35.0±3.0 


9.7 ± 1.6 


LM88 



[b] 

Table 2. List of 37 pulsrs which are analysed in Figure 2 and 3. 
The corresponding emission altitudes presented in Fig. 2 and 3 
are taken from Paper I and II, as well as from the Table 1 in this 
paper. 



B0329±54 


B0355±54 


B0450-18 


B0525±21 


B0540±23 


B0740-28 


B0823±26 


B0833-45 


B0834±06 


B0919±06 


J1022±1001 


B1133+16 


B1237±25 


B1642-03 


B1706-16 


J1713±0749 


B1727-33 


B1749-28 


B1800-21 


B1804-08 


B1823-13 


B1845-01 


B1848±04 


B1855±09 


B1859±03 


B1913+16 


B1915+13 


B1916+14 


B1924±16 


B1933±16 


B1952±29 


B2020±28 


B2021±51 


B2310±42 


B2319±60 


B2324±60 


B2334±61 









ation at a given frequency v is generated) can be expressed 
in the form 




where p is the opening angle [in degrees] between the mag- 
netic axis and the tangent to dipolar magnetic field lines 
at points where the emission corresponding to the appar- 
ent pulse- width originates, and s is the mapping parameter 
(0 < s < 1) describing the locus of field lines on the po- 
lar cap (s = at the pole and s = 1 at the edge of the 
Goldrech- Julian (1969) circular polar cap). Opening angles 
p — p(W;a,/3) are calculated using the knowledge on the 
viewing geometry (represented by the inclination angle a 
between the rotation and the magnetic axes and the im- 
pact angle /3 of the closest approach of the observer to the 
magnetic axis; see Paper I and references therein for more 
detailed explanation), as well as measurements of the pulse 
width W. Estimates of emission altitude (for s « 1) from 
equation (2) are listed in Table 1. 

In Figure 1, we present a plot of emission altitude r cm 
versus period P for very young and very old pulsars cor- 
responding to pulse width measurements at 1.4 GHz. The 
data have been taken from Table 1 and from the Paper II (for 
MSPs: J1022+1001, J1713+0747 and B1855±09). It is no- 
ticeable that emission altitudes differ slightly between very 
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young and very old pulsars. A formal weighted fits to very 
young and very old pulsars yield r em = i?(55±8)P~°' 25±0 ' 08 
and r cm = R(35 ± 6)p-°- 38±0 04 , respectively. We can now 
gather all estimates of emission altitudes from Papers I, II 
and from this paper (Table 1) to obtain a plot of emission al- 
titude versus characteristic age for 37 pulsars (see Fig. 2 and 
Table 2) . The periods of normal pulsars were divided into six 
groups: (0.05 s- 0.2 s), (0.2 s - 0.4 s), (0.4 s - 0.6 s), (0.6 s 
- 0.8 s), (0.8 s - 2.0 s), (2.0 s - 4.0 s) and data belong- 
ing to each group were represented by different size of the 
circle (with biggest circles corresponding to longest periods 
in our sample). This is consistent with equation (1), show- 
ing that for a given characteristic age, pulsars with longer 
periods emit their radio emission from higher altitudes than 
those with shorter periods. A formal weighted fit to all data 
points in Fig. 2 is r om = i?(33 ± 1)t^~°' 07±0 ' 02 , confirming 
a putative weak dependence r cm on pulsar age r = re • 10 6 
years (equation 1). 



3 DISCUSSION 

In Papers I and II the radius-to-frequency mapping and the 
period dependence of emission altitude have been discussed 
in details. Here, a detailed analysis of age dependence on 
emission altitude is carried out. Results presented in Figs. 1 
and 2 indicate that emisson altitudes depend on the charac- 
teristic age r C m(ra) tx t^ ' 07 , which is consistent with equa- 
tion (1). In this analysis we used 37 pulsars at 1.4 GHz for 
which all necessary information were available, that is: very 
good quality profiles (pulse width measurements) as well as 
reliable estimates of a and /3 angles. 

It is important to notice that the radio emission re- 
gion for a group of very old pulsars (containing millisecond 
and normal pulsars) is located at correspondingly lower al- 
titudes than that of very young objects (see Figs. 1 and 2). 
Three points in Fig. 2 deviate significantly from the rest of 
distribution. These points correspond to 3 millisecond pul- 
sars in our sample. One could think that this is an effect 
of nondipolar magnetic field in the radio emission region of 
millisecond pulsars. However, we do not regard these outlier 
points as strong indication of deviation from dipolar fields. 
In fact, the emission altitude at a given frequency is deter- 
mined more strongly by the value of pulsar period than by 
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Figure 2. Emission altitudes versus characteristic age r for 37 pulsars. The size of circles represents increasing period from 0.1 s to 4.0 : 
(see text for explanation). The filled dots correspond to millisecond pulsars. A formal weighted fit gives r om = R(33 ± l)T e r ' 07±0 ' 02 . 



its characteristic age (see equation 1). Therefore, the emis- 
son altitude of these milisecond pulsars is about one order 
of magnitude smaller than in typical pulsars, with periods 
in the range 0.1 — 4.0 seconds (thus covering one order of 
magnitude), much longer than the millisecond ones. If one 
omits three points corresponding to the millisecond pulsars, 
then a formal weighted fit is r cm = i?(35±l)T^~ 04±02 , again 
consistent with equation (1). 

The surface magnetic field is estimated using the as- 
sumption that the main rotational energy loss is in the form 
of electromagnetic dipole radiation. With a few further as- 
sumptions about the pulsar spindown mechanism (Manch- 
ester & Taylor 1977; Camilo, Thorsett & Kulkarni 1994), 
the surface magnetic dipole field component of a pulsar (in- 
ferred from its values of P and P) can be expressed in terms 
of pulsar characteristic age r = 0.5 • P/P [s] as follows 



B d =4x io 12 -rr 1/2 



P [G], 



(3) 



where ra 



s/10 6 . If the magnetic field in the emission 



region is dipolar, one can calculate its values at the emission 
altitude r cm as 



B T — B, 



( rem V 
V R ) 



From equations (1), (3) and (4), one obtains 
Bl « 2 x 10 7 4 H 3 z r 6 - - 29 [G]. 



(4) 



(5) 



The above equation indicates that the magnetic field in the 
emission region may depend on the characteristic pulsar age. 
In typical pulsars with re ~ 1, the magnetic field in the 
emission region corresponding to 1.4 GHz should be about 
10 7 G. 

Figure 3 presents a plot of emission altitude versus the 
dipolar component of surface magnetic field (calculated from 
equation 3) and the magnetic field in the emission region 
(calculated from equation 4). The magnetic field in the emis- 
sion region at 1.4 GHz for normal pulsars (10 6 — 10 7 yr) has 
values of about 10 7 - 10 s G. Gil & Kijak (1993) suggested 
that the magnetic field in the emission region B v is almost in- 
dependent of the pulsar period and characteristic age. This 
analysis clearly shows that very old pulsars (filled circles) 
and very young pulsars (circles with stars) have quite dif- 
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Figure 3. Plot of the emission altitude r em versus two kinds of magnetic field, the surface magnetic dipolar field component B<j (squares) 
and the field strength at the estimated emission altitude B r (circles). The filled circles correspond to very old pulsars. Circles with star 
correspond to very young pulsars. The PSRs B1534+12 and B1913+16 are marked by the circles with horizontal arrows (see Section 3). 



ferent magnetic field in the emission region B T (Fig. 3 and 
equation 5). 

Recently, a relationship between the characteristic and 
true age was discussed in the literature. Several authors 
(Lyne, Anderson & Salter 1982; Camilo, Thorsett & Kulka- 
rni 1994; Bhattacharya 1996; Lorimer 1996, Cordes & 
Chernoff 1998) suggest that the true age is somewhat 
smaller than the characteristic age computed from P and P 
measurements. Arzoumanian, Cordes & Wasserman (1999, 
hereafter ACW) argued that the three typical pulsars 
PSRs B1534+12, B1913+16 and B2127+11C are probable 
younger than their characteristic age. Let us consider two 
of them, PSRs B1534+12 and B1913+16, for which the val- 
ues of a and f) angles are known (Arzoumanian et al. 1996; 
R93, respectively). In Table 3, the geometrical emission alti- 
tudes r e m(W;a,/3), as well as semi-empirical altitudes tkg 
calculated from equation (1) for characteristic age (for P, 
P) and "true age" (taken from ACW), respectively, are 
listed. One can see that rKGrj6(ACW)] are closer to r em 
than rKG[Te{P, P)]- Moreover, as one can notice from Fig. 



3, these two pulsars seem to have B r values considerable 
lower than the rest of the group of normal pulsar (between 
two dashed vertical lines) . It is worthy of noticing, that even 
their characteristic age (not to mention the "true age") is 
much smaller than rest of the old pulsars (filled circles). 
However, if we use "true age" (after ACW) in equation (5) 
then the verified values of B T increase such that they fall 
into a region between dashed vertical lines (which is marked 
by horizontal arrows). This indicates that the emission alti- 
tudes r om and hence the magnetic field B T in the emission 
region depends on the true age rather than on the charac- 
teristic age. Of course, in most cases, when the actual period 
P is much larger than the initial period, these values are the 
same. 

Finally, let us discuss a possibility that the magnetic 
field lines in the emission region are not purely dipolar. The 
observational opening angle p — p(W;a,/3) is calculated us- 
ing geometrical equation proposed by Gil, Gronkowski & 
Rudnicki (1984), which holds on the pulsar centered ce- 
lestial sphere, independently of geometry of magnetic field 
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Table 3. Emission altitudes for two PSRs: B1534+12 and B1913+16 at 1.4 GHz. Data for the geometrical analysis arc taken from 
Paper I (B1913+16) and Arzoumanian, Phillips, Taylor & Wolszczan (1996) for B1534+12. Emission altitudes of tkg are calculated for 
characteristic age and age proposed by Arzoumanian ct al. (1999). 



PSR 


Geometrical 


Characteristic age 


"True age" 




analysis 


analysis 


analysis 


B1534+12 


r om = 22.R 


t 6 = 252, r KG (252) = 12R 


r 6 < 10, r KG (10) > 15fl 


B1913+16 


r cm = 18fl 


T(j = 109, r KG (109) = UR 


t 6 < 5, r KG (5) > 18R 



lines in the emission region. This observationally inferred 
opening angle p — p b B {W; a, f3) is also the angle between 
magnetic axis and the tangent to magnetic field lines at 
points where the emission correspondig to the apparent 
pulse width W originates. For dipolar field lines p is de- 
scribed by equation (2), with s — 1 for last open field lines 
corresponding to low level emission below 1% of the max- 
imum intensity (e.g. Table 1). For a simple model of star- 
centered quadrupole field in the radio emission region, the 
opening angles p q = 0.8s(r q / R)P~ X ^ 2 deg (where s < 1), 
as compared with p d = 1.2(rd/-R) 1//2 -P~ 1//2 deg for dipolar 
geometry. Assuming that p d = p q = p obB (W; a, f3) we get 
ry 2 /r q = s/R 1 ^ 2 , where and r q represent the emission 
altitudes for dipole and quadrupole field lines, respectively. 
Thus, if one allows nondipolar magnetic field in the emission 
region, then equation (1) represents the upper limit for emis- 
sion altitudes. Moreover, the radius-to-frequency mapping 
would change from r oc v~ v for dipolar field toroc v^ v ^ 2 for 
quadrupole field, where the index p is constrained observa- 
tionaly as 0.2 ±0.1 (see Table 4 in Paper II). Unfortunately, 
this constraint is too weak to differentiate between dipole 
and quadrupole field lines. 

Theoretically, a deviations from dipolar form can be 
caused by (i) outflow of relativistic currents, (ii) field lines 
sweep-back due to fast rotation, and (iii) currents flowing in 
the thin layer of the neutron stars outer crust. Let us briefly 
discuss these three cases: 

(i) The outflow of the charged component of the magneto- 
spheric plasma cannot change a structure of strong magnetic 
field generated inside the neutron star. In fact, using the 
Ampere's law § B t ■ ds = (4tt/c)J p , where B t is the toroidal 
component generated by the maximum available poloidal 
current through the polar cap J p = 47rcp G .jr 2 (r cm /i?) 3 ^ 2 , 
where r cm > R is the emission altitude, r p « 10 4 .P _1,/2 cm 
is the polar cap radius, p G j = B/cP is the Goldreich- 
Julian (1969) corotational charge density, c is the speed 
of light, B is a poloidal (dipolar) pulsar magnetic field 
and P is the pulsar period. One can easily show that 
B t /B « 10" 6 P~ 3/2 (r cm /R) 3/2 . Thus, for typical pulsar B t 
is absolutely negligible, while for the millisecond pulsars the 
toroidal contribution to dipolar magnetic field can achieve 
a few percent at most. 

(ii) The sweep-back of field lines by corotation with neu- 
tron star can occur near the light cylinder, where the coro- 
tation velocity approaches the speed of light. However, it is 
generally accepted that the radio emission originates at al- 
titudes much smaller than the light cylinder radius (Cordes 
1978, BCW, Kramer et al. 1997). 

(iii) The currents in the crust are most probably sources of 



multipolar components of the pulsar magnetic field outside 
the neutron star. However, these contributions decay rapidly 
with increasing distance from the star surface, and at radio 
emission altitudes only dipolar field can be expected, at least 
in normal pulsars. In the millisecond pulsars, in which the 
emission region can be as close as few stellar radii from the 
polar cap, the magnetic field lines can, in principle, devi- 
ate significantly from dipolar form. However, Arons (1993) 
has shown that the location of the spin-up line in the P — P 
diagram constrains the surface strength of multipole compo- 
nents to no more than 40% of the dipole field in millisecond 
pulsars. 



4 CONCLUSIONS 

The main conclusions from this paper and from previous 
work (Papers I and II) on emission regions in radio pulsars 
are following: 

(1) A radius-to-frequency mapping operates in the emis- 
sion region and has a form of r cnl (u) oc j,- () - 21 ± - 07 

(2) Pulsar radio emission is typically generated at alti- 
tudes smaller than a few per cent of the light-cylinder ra- 
dius tlc and the ratio t*k G /?"lc decreases with increasing 
period. The total extent of the emission region is smaller 
than 250.R in longer-period pulsars and correspondingly less 
in shorter-period pulsars. 

(3) The emission region in old pulsars is located at cor- 
respondingly lower altitudes than in young pulsars with ap- 
proximately the same period. 

In summary, pulsar emission altitudes depend on fre- 
quency v, pulsar period P and pulsar age r, as described by 
equation (1). 
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